Effect of electron correlations in Pd, Ni, and Co monowires 
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We investigated the effect of mean-field electron correlations on the band electronic structure 
of Co, Ni, and Pd ultra-thin monatomic nanowires, at the breaking point, by means of density- 
functional calculations in the self-interaction corrected LDA approach (LDA-I-SIC) and alternatively 
by the LDA-l-17 scheme. We find that adding static electron correlations increases the magnetic mo- 
ment in Pd monowires, but has negligible effect on the magnetic moment in Co and Ni. Furthermore, 
the number of d-dominated conductance channels decreases somewhat compared to the LDA value, 
but the number of s-dominated channels is unaffected, and remains equal to one per spin. 

PACS numbers: 75.75.-(-a, 73.63.Nb, 73.22.-f, 75.20.En, 75. 15. Mb 



I. INTRODUCTION 

Nano- and molecular-sized components hold consider- 
able promise as forming the platform for an entirely new 
type of hardware technology. In such a future technology, 
quantum effects inherent in these nano-objects might give 
rise to new functionalities, possibly resulting in a totally 
different design of the basic components in information 
technology from what we are used to today. An impor- 
tant category of systems is constituted by ultra-thin con- 
ductors and nanocontacts. The best known and most 
investigated systems in this class are carbon nanotubes 
and molecular conductors. Metallic nanowires are also 
of high interest, but much less studied. Ultimately thin 
nanowires, i.e. freestanding atomic chains hanging be- 
tween tips, constitute a unique metallic structure where 
nanosize properties of matter can be tested. Unlike nan- 
otubes, they are clearly transient objects, although for 
applications they could be stabilized, e.g. inside other 
nanostructures. Besides, their transient nature could it- 
self be put to use. 

Metallic nanowires are fabricated using a scanning tun- 
neling microscope (STM) or an atomic force microscope 
tip. By forming a metallic contact between the tip and 
sample of the same metal as the tip, and then retract- 
ing the tip, a thin neck or wire may form by plastic de- 
formation. Just before it breaks, the nanocontact may 
consist of a single atom, and under some circumstances 
and in some metals of a short single chain of atoms - 
a segment of monatomic nanowirei^ A second way in 
which nanowires can be produced is inside a transmission 
electron microscope- (TEM). With the TEM beam, two 
nearby holes are burned in a thin metal film, close enough 
that the bridge between them narrows down sponta- 
neously to form a nanowire. A third way to produce 
nanowires is via the mechanically controllable break junc- 
tion techniquei^ 

It was recently predicted that quantum confinement 
of the electronic states in metallic atomic chains should 
lead to a magnetic ground state, and to a superpara- 



magnetic state at finite temperatures, in Ru, Rh, Pd, 
Os and ptiSiLLS In the present paper, we build on 
those results, also stimulated by some recent intriguing 
experimental results indicating unexpectedly low frac- 
tional conductance through Co, Pd and Pt nanocon- 
tacts and nanowires observed at room temperature^iifliii 
(note however that at low temperature, the uncontami- 
nated break junctions fail to show fractional conductance 
and even nanowires in the case group III and IV transi- 
tion metals.^") Assuming that, as in the room temper- 
ature experiments, short nanowires could easily form at 
the nanocontact, we wish to explore the possible effects 
of strong correlations on the electronic structure of the 
wires, and if and in what way strong correlations might 
relate to these recent experimental results. 

We will be concerned here only with those correla- 
tion effects that can be included within conventional 
Fermi liquid theory, where quasiparticle electronic bands 
replace the noninteracting electron bands, but can be 
treated as such in every other aspect. Thus, we will not 
be addressing at all the strongest kind of correlation ef- 
fects, including Coulomb blockades and Kondo phenom- 
ena, probably more typical for quantum dots than for 
genuine metallic nanocontacts. The effects which we will 
address here still amount to a nontrivial renormalization 
of the electronic bands, which can however be treated by 
means of static mean-field theories such as LDA-I-C/ and 
LDA+SIC. 

The electron bands in transition metal nanowires, es- 
pecially the d bands, are considerably narrowed by the 
reduced one-dimensional coordination. Band narrowing 
is further enhanced under the strong tensile stress that 
brings the nanowires close to the point of breaking. It is 
well known that for narrow band systems, standard den- 
sity functional theory (DFT) approaches'"^ - notably the 
local density approximation (LDA) and generalized gra- 
dient approximations (GGA) - do not reproduce well the 
electronic localization within the polarized atomic shells. 
For that reason we find it relevant and instructive to 
investigate how extensions of DFT, devised specifically 
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to improve the mean field static description of electron- 
electron correlations, affect the electronic structure of 
strongly stretched monatomic wires. In particular, we 
will make use of the LDA-I-J7 approachi^ and the self- 
interaction corrected LDA (LDA-f-SIC) scheme.^"* 

We may anticipate that with these corrections, the 
band structure near the Fermi level Ep, including band- 
widths and splittings could be strongly altered. In that 
case also the channel number (number of propagating 
bands that cross Ep), their transmission, and in the end 
ballistic conductance might be sensitive to correlations. 
In this work, devoted to strictly ideal infinite monatomic 
nanowires, we will address the sensitivity of bands and 
channels, though clearly not of transmission, to strong 
(albeit for the time being statically described) electron- 
electron correlations. 



II. METHOD 

The geometry adopted for the calculation was a two- 
dimensional lattice of infinite ideal monatomic chains, 
with the spacing between chains large enough to make 
their interaction negligible. 

We performed all calculations using the SIESTA- 1.3 
codeii^ The starting point is in all cases a local (spin) 
density approximation, L(S)DA, calculation. On top 
of that, we treat the effect of strong correlations by 
an extension of LDA, namely the LDA+U methodic 
and the self-interaction (SI) corrected LDA approach 
(LDA-f SIC).^"' The latter method is an approximated 
scheme to retain the exact exchange and it does not add 
any correlations; LDA-I-SIC even removes a part of cor- 
relations, i.e. self-correlations. Thus, the effect of SI 
correction may differ from the LDA-f-U result. 

The LDA-I-U correction potential, VjJ^^, to the LDA 
Hamiltonian is diagonal, and defined as 

(u-j) J2 «' - o + 

{U-J){l-n-,), (1) 

with being the average orbital occupation of the cor- 
related shell calculated self-consistently from the orbital 
occupations n^. 



from the MuUiken occupations n^^""**^'^" (Ref. [p) as 
follows 
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during the LDA+[/ iteration process. We use here for Nj, 
the total actual occupation number of the polarization 
shell (for instance d-shell) rather than fixing nJJ at half 
occupation or using an average LDA occupation. Details 
and a test of this implementation are given in Ref. |l(il |. 
Our definition of the local occupations also differs 
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involving differently the density matrix Df^i, and the over- 
lap integrals Smfi = ("il/x) and S^m = Wl'm) for the basis 
functions /i and u and the projector function m. Further 
discussion of the local occupation numbers is given in 
Appendix A. 

The d-shell exchange parameter J was set to 1 eV in all 
calculations. Varying J in the range 0.5 — 1.5 eV was seen 
to have only a very minor effect on the studied systems. 
Several values for the intra-atomic Coulomb parameter U 
in the range 3 — 9 eV were tested in order to get a broad 
insight into the electron localization effects. Thus, we 
use J7 as a free parameter in our calculations, although 
it can in principle be calculated self-consistently. We 
did calculate U self-consistently for our nanowires using 
the constrained density functional method described in 
detail in Refs. and 0, in addition to treating U 
as a free parameter. The results are detailed and crit- 
ically assessed in Appendix A. However, an indiscrim- 
inate use of self-consistently calculated U parameters is 
not unproblematic. The LDA-f [/ method has indeed had 
great success in describing systems with rather localized 
d electrons, for example the 3d oxide NiO, where the 
LDA-f [/ gives the correct insulating ground state as op- 
posed to standard LSDA. On the other hand, when the 
LDA-I-J7 method is applied to systems which are in ef- 
fect not strongly correlated, for example TiO, it fails, at 
least when a self-consistently calculated U is used in the 
calculation, 

In the present implementation of the LDA+SIC 
scheme we adopted a modified pseudo-SIC approach 
closely akin but not identical to that in Ref. [l^ . 

The LDA+SIC correction matrix, Vgj(j to the LDA 
hamiltonian is given by 



cr 



J2 (Ml7f)(7f 1^) 



= V, 



Cf 

K(r)]0.(r), 



\V^xcmr). 



HXC 



(5) 

(6) 
(7) 



with pf (r) being the orbital density and the projectors 
(/)i(r) being the same basis functions /i, v used for the 
calculation of local occupation numbers in the LDA+J7 
method. Above Sl-correction potential differs from the 
original potential introduced by Filippetti and Spaldin by 
a factor of 1/2 which we dropped in order to obtain the 
exact result in the one-electron limit {e.g. the Hydrogen 
atom). The second difference between our approach and 
the LDA-hSIC in Ref. [13 is the fact that, we use as a 
projector the shortest basis function (last zeta-function) 
instead of long-range pseudopotential. In present work. 
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we restricted the self-interaction correction to the d-shell 
only. 

SIESTA treats core electrons by employing the 
norm-conserving pseudopotentials in their fully non- 
local (Kleinman-Bylander) formi2£ We generated the 
relativistic pseudopotentials for atomic elements us- 
ing the Martins- Troullier— scheme with the nonlin- 
ear core corrections.— The electronic configurations 
and the cut-off radii (in a.u.), in the s/p/d/f or- 
der, were: As^Ap^Sd'^if and 2.0/2.0/2.0/2.0 for 
Co, 4524^03^84/0 and 1.7/1.88/1.88/1.88 for Ni, 
5515/4^94/0 and 2.0/2.2/2.0/2.2 for Pd in the wire and 
2.39/2.39/1.79/1.79 for Pd in the bulk. 

The basis set adopted is a very general and flexible lin- 
ear combination of numerical atomic orbitals (LCAO)^ 
by Sankey and NiklewskL^"* It allows for arbitrary angu- 
lar momenta, multiple-C, polarized and off-site orbitals. 
For each atom, we used an automatically generated dou- 
ble zeta basis set with polarization functions (which are 
p-type functions in the case of transition metals). The 
second zeta of the d-shell was used as projector function 
for the LDA+U and LDA-I-SIC operators. The radii of 
these projectors were (in a.u.): 1.95 for Ni, 2.073 for Co, 
and 2.675 for Pd. For the self-consistent calculation of 
the parameter U we tested in addition projector radii of 
1.8 a.u. for Ni, and 2.2 a.u. for Pd. 

As for the Monkhorst-Pack gridf^ for the bulk calcu- 
lations we used 20x20x20 mesh-points, and for all wires 
one k-point in the plane perpendicular to the wire and 
100 k-points along the wire. This is good enough to rep- 
resent a single nanowire in this case, in case of a very 
large separation between nearest neighbour nanowires. 
Accordingly, the distance between wires in the perpen- 
dicular plane was set to 10 A. For the real space grid, we 
used a uniform mesh of quality corresponding to an en- 
ergy cut-off of 400 Ry in the plane- wave method (which, 
in case of wires with one atom per cell, results to 2187000 
mesh-points) . 



III. RESULTS AND DISCUSSION 

A. Interatomic distances and magnetization within 
LSDA 

The top panel of Figure ^ shows the LSDA total en- 
ergy as a function of interatomic distance. Our SIESTA 
calculations are denoted by closed symbols and (1) in the 
legend. We see that for all three metals, the total energy 
has a well-defined minimum, which we will call the equi- 
librium interatomic distance. For both Ni and Co, this 
distance is around 2.1 A, whereas for Pd, it is significantly 
larger, around 2.5 A. Decreased coordination from 12 to 
2 is of course responsible for the considerable decrease 
from bulk interatomic distances, 2.42 , 2.42 and 2.56 A 
respectively for Ni, Co, and Pd. For comparison, we also 
show total energy curves obtained by other electronic- 
structure codes (using LDA), denoted by open symbols 



and (2) in the legend. For Co and Ni we plotted in par- 
ticular results obtained by Smogunov et al. in Ref. 
by means of plane wave method using the PWscf code^ 
for Pd we show results of an all-electron FP-LMTO (see 
Ref. 26] for the code) calculation, independently per- 
formed by us. As is seen, the agreement between results 
using the SIESTA code and other codes is more than 
reasonable. 

The ideal breaking point a;, of a nanowire is defined 
as the atomic distance at which a wire under increasing 
stress ceases to remain locally stable and must disappear 
(if it did not break earlier at the wire-tip junctions as 
it happens in practice). This is an important parameter 
since in a tip-wire-tip nanocontact the very last conduc- 
tance plateau before the wire breaks under increasing 
strain will refer to interatomic distances that are larger 
than equilibrium due to stress, but still have as an 
upper limit. 

As defined, the ideal breaking point a;, will coincide 
with the smallest interatomic distance associated with 
instability of some normal mode, typically with a phonon 
softening. While calculating the full phonon spectrum of 
the wire is beyond the scope of the present work, we can 
crudely estimate ah in a rather simple way. In a chain 
held together by pairwise short range forces the longi- 
tudinal vibration branch is unstable when for increasing 
spacing a the interatomic restoring force is maximally 
negative, or equivalently, when the second derivative of 
the total energy is zero. Although metal nanowires are 
not exactly a case of pairwise short range forces, the 
crude estimate obtained in that approximation is suffi- 
cient for our purposes. In the middle panel of Figure ^ 
we plot the tension force, F, defined as 

F = pA, (8) 

where p is pressure per cell and A is the surface of a plane 
perpendicular to the wire and restricted by the walls of 
the unit cell (in our case this is a hexagon with an area 
of 541.265877 A^, so the pressure of 1 kBar in a cell with 
1 atom corresponds to the tension force of 5.68 nN). The 
force in the Pd wire is seen to have a local minimum 
around 2.9 A, and the corresponding values for Ni and 
Co are both 2.60 A. In the following, we will identify 
these distances as the ideal breaking points of the wires, 
keeping in mind that as explained above, they are re- 
ally extreme and rough upper limits of the true breaking 
points. We note also that for a Pd monowire the force 
curve is very flat, indicating that the interatomic spacing 
in that system under stress should be very fluctuating. 

The bottom panel of Figure ^ shows the spin mag- 
netic moment per atom in //^ as a function of the inter- 
atomic distance. For Pd, we also plot the magnetization 
calculated using the FP-LMTO code£ The magnetiza- 
tions obtained for Ni and Co wires (flUed symbols) are 
compared with the PWscf calculations given in Ref. 
(open symbols). The magnetization results are summa- 
rized in Table ^ where we compare them with bulk in- 
teratomic distances and magnetic moments. 
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TABLE I: Bulk and nanowire interatomic distances and 
magnetic moments calculated with SIESTA in LSDA. The 
bulk calculations were performed assuming an ideal hexago- 
nal structure for Co, and fee for Ni and Pd. 



Interatomic distance [A] 



FIG. 1: Total energy, tension force and magnetization per 
atom for Co, Ni and Pd monowires as a function of in- 
teratomic distance. (1) denotes results obtained with the 
SIESTA method and (2) denotes results obtained with other 
methods; references are given in the text. The total energy 
curves for Co, Ni and Pd have been shifted, so that the min- 
imum is at zero for all curves. 



The magnetic moments for Ni and Co calculated with 
SIESTA are somewhat higher than these numbers calcu- 
lated with the plane wave method, and than experimental 
bulk magnetizations. The onset of magnetism for Ni and 
Co wires occurs for slightly smaller interatomic distances 
in the SIESTA calculations compared to the PWscf calcu- 
lations. For the Pd wire the magnetism occurs at slightly 
larger interatomic distances and appears to be more grad- 
ual than the magnetic moment calculated with the FP- 
LMTO code. The source of these slight discrepancies 
is related to differences of used pseudopotentials and the 
difference between the PW and LMTO schemes. SIESTA 
uses the normconserving PPs unlike the PWscf calcula- 
tions where ultrasoft PPs were employed. We checked 
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Co 


Ni 


Pd 


interatomic distance (A) 






bulk, this work 


2.42 


2.42 


2.56 


bulk, exp.29 


2.51 


2.46 


2.67 


wire equilibrium distance 


2.08 


2.08 


2.50 


wire breaking distance, 


2.60 


2.60 


2.90 



magnetic moment (/is) 





in bulk, this work 


1.95 


0.63 


0.0 




in bulk, expf22i 


1.72 


0.61 


0.0 




at wire equilibrium distance 


2.23 


1.21 


0.45 


J 


at wire breaking distance 


2.46 


1.32 


0.69 


- 


free atom configuration 









also the effect of using the pseudopotentials with and 
without nonlinear core-corrections^^ (NLCC) for mag- 
netism. For the PPs with NLCC included, magnetism 
in Pd wire occured at the interactomic distance larger of 
about 0.2 A. 

All above calculations were performed with unit cells 
containing only one atom, therefore the ferrmomagnetic 
order in wires was assumed. 

Finally, we tested the possibility of antiferromagnetic 
ordering in Pd, Co and Ni wires. Antiferromagnetic (AF) 
exchange could in principle prevail over direct ferromag- 
netic (FM) exchange in the extreme tight binding limit 
attained for very large a. We performed the AF and FM 
calculations with two atoms in the unit cell (to compare 
the total energy of two phases). We also calculated the 
chain of seven atoms in the unit cell, flipping the spin 
in the middle atom. We found that, coUinear antiferro- 
magnetic ordering is generally unstable with respect to 
ferromagnetic ordering, at the breaking point distance as 
well as for shorter interatomic distances, and in all metals 
studied in this work. 



B. Conductance channels and magnetic moments - 
effect of electron correlations 

A central aim of the present work is to elucidate the 
effect of correlations on the number of conductance chan- 
nels, foreshadowing in a qualitative way the effect on bal- 
listic conductance. In the Landauer picture of nanowire 
conductance, each nanowire electron band crossing the 
Fermi level constitutes a conducting channel between 
the tips at the two ends of the nanowire. Each channel 
can transmit electrons fully or partially, depending on a 
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FIG. 2: Nanowire band structure with orbital labels (upper 
panel) and DOS (lower panel) for the majority spin of Pd 
obtained with LSDA at the breaking point. The Fermi level 
is at zero. The symmetry point A in the band-structure plot 
corresponds to the point (0,0, | ^) in the Brillouin zone. 



transmission coefficient dictated by the tip-wire junction, 
and by the character of the electron band. In transition 
metals, bands with s character are generally close to total 
transmission. Owing to their relatively large bandwidth, 
s electrons are little disturbed by the potential change 
they experience at the tip-wire junction, and are thus 
essentially not reflected back, with transmission above 
95%. In d-dominated channels, on the other hand, the 
smaller bandwidth of d states leads to higher reflection at 
the junction, with a much lower transmission, typically 
around 20% (Ref. Q). 

The d bands should of course be narrowest, and corre- 
lations strongest, for a nanowire under maximal strain. 
Thus, in order to investigate the effect of electron cor- 
relations on the number of conductance channels of the 
wires, we focus here on the nanowires electronic band 
structures at the breaking points. These increased inter- 
actomic distances are more similar to the conditions at 
which the experiments on breaking junctions have been 
done. Further, we identify the amount of s and d charac- 
ter of the bands at the vicinity of the Fermi level. This 
enables us to draw some qualitative conclusions regarding 
the conductance of the wires. In particular, it will allow 
a discussion of whether or not a conductance significantly 
below the unit conductance quantum Go = 2e^//i is or is 
not directly suggested by electron correlations, at least 



within the present static description. 

We start by identifying the character of the bands 
around the Fermi level by comparing the majority spin 
orbital-character resolved partial density of states (DOS) 
with the band structure. In Figure |21 we show the band 
structure and s, d^^, dxy + dj.2_y2, and dxz + dyz par- 
tial DOS for the Pd wire at the breaking point, obtained 
from the LSDA calculations. In the band-structure plot, 
the character of the individual bands has been indicated 
based on information in the DOS plot. When the band 
is s — c? hybridized but one of the contributions is much 
smaller than the other, we show the smaller component 
within brackets. The band structures for Co and Ni are 
generally similar to that of Pd, and we refrain from re- 
peating the same analysis over and over. For Co and 
Ni the character of each band can be obtained by direct 
comparison with Pd. 

Figure 121 shows the electronic band structures for the 
Co, Ni and Pd wires at their respective breaking points. 
The left-most column shows first the standard LSDA cal- 
culation, the four middle columns the LDA-|-[/ results for 
a range of U between 3 eV and 9 eV, and the right-most 
column the results obtained in LDA-I-SIC. The overall 
effect of added correlations on the electron structures is 
that bands shift in energy and that the dispersion of the 
s + d^i bands change in the vicinity of the Fermi level. 
The forms of the individual bands, however, are to a large 
degree conserved. We also see that the effects of includ- 
ing self-interaction corrections on the band structure are 
roughly similar to the effects of the LDA+[/ calculation 
with U around 3 eV or smaller. 

There are quite substantial differences in the way 
LDA-I-C/ and LDA+SIC affect the electronic structure. 
The SIC potential is linearly dependent on occupation 
numbers. Therefore unoccupied bands do not shift, and 
nearly empty bands shift very little, as for example the 
minority spin d^y + dx2-y2 band in Co. On the con- 
trary, the LDA-I-J7 approach moves unoccupied bands 
upwards with a constant shift {U — J)/2, due to the last 
term in Eq. , (see Section llT|) . As we mentioned in the 
SectionHU the LDA+SIC and LDA+U results may dif- 
fer. This is due to the fact that LDA-I-SIC does not 
add correlations; the method removes "self-exchange" 
and "self-correlations". While the LDA-I-C/ method is 
self-interaction free within the considered shell and addi- 
tionally adds correlations within this shell. 

The magnetic moments of Co and Ni are largely unaf- 
fected by increased correlations, whereas the Pd moment 
increases substantially with U, see Table Ull 

In general, the properties of transition metal systems 
are strongly sensitive to the relative position of the s and 
d bands and the charge transfer between these two bands. 
For example, the magnetism in Pd wires is governed by 
d ^ s charge transfer.^ Therefore, it is of interest to 
study the effect of correlation on the d s transfer. 
Table ^ shows that the d — > s charge transfer decreases 
with U for Pd (and thus in principle works toward killing 
off the Pd magnetism by opening up a gap), but it is 
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FIG. 3: Band structures of Co, Ni and Pd nanowires obtained in LSDA; in LDA+U for U = 3.0, 5.0, 7.0 and 9.0 eV and J=1.0 
eV; and in LDA+SIC (d-shell only). Thin line: spin up electrons; thick line: spin down electrons; dashed line: Fermi level. All 
bands are plotted at the breaking point interatomic spacing ai,. 



TABLE II: Properties of monatomic wires at the breaking point. AIt is the magnetic moment per atom in hb, Ma the 
spin-polarization of the s shell in /is, nj"*"-^ the total MuUiken occupation of the s shell, and Nc the number of spin-resolved 
conductance channels. 



Method 


Cobalt 




Nickel 


Palladium 






Mr Ms 




Mt Ms nl+^ (Nl,Ni) 


Mt Ms 





LSDA 




2.46 


0.03 


1.40 


(1,6) 


1.32 


0.03 


1.24 


(1,6) 


0.69 


0.29 


0.50 


(1,4) 


LDA+U (3 


eV) 


2.44 


-0.01 


1.43 


(1,6) 


1.32 


0.03 


1.24 


(1,6) 


0.69 


0.29 


0.50 


(1,1) 


LDA-HU (5 


eV) 


2.43 


-0.04 


1.48 


(1,3) 


1.37 


0.08 


1.20 


(1,4) 


0.84 


0.30 


0.62 


(1,1) 
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2.44 


-0.04 


1.48 


(1,3) 


1.42 


0.07 


1.23 


(1,4) 


0.96 


0.22 


0.80 


(1,1) 


LDA+V (9 


eV) 


2.45 


-0.03 


1.47 


(1,3) 


1.44 


0.05 


1.26 


(1,4) 


0.97 


0.13 


0.89 


(1,1) 


LDA-hSIC 




2.43 


0.01 


1.42 


(1,6) 


1.28 


0.02 


1.23 


(1,5) 


0.56 


0.24 


0.40 


(1,1) 



much smaller for Ni and roughly constant for Co. 

In the Co wire, the largest change of the band struc- 
ture when correlations are added is the upward displace- 
ment in energy of the flat nondispersive minority spin 
dxy + (1^2 _y2 band. The minority channel of this orbital 
sits right at the Fermi level in the LSDA and in LDA-I-SIC 
calculations, and it is pushed up from the Fermi level with 
increasing correlations within the d-shell in LDA-|-?7. For 
U=7 eV and above, it is no longer visible in the band- 



structure plots. In contrast to that, the majority spin 
channel of this orbital is pushed down in energy. All re- 
maining minority spin orbitals shift slightly down, so as 
to accommodate the charge transfer from the emptied 
minority -|- (ia;2_j^2 orbital. The corresponding major- 
ity spin bands are also downward shifted. For U—3 eV, 
at the band edge at A, the d^^ + (s) band of the minority 
spin channel has moved down slightly and positioned it- 
self just on the Fermi level. For C/=5 eV and larger, four 
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FIG. 4: Spin- and orbital-resolved s and d^2 DOS for the Pd 
wire at the breaking point, obtained with LSDA (top panel), 
LDA-I-C/ (four middle panels) and LDA+SIC (bottom panel). 
The Fermi level is at zero. 



bands cross the Fermi level: the minority spin d^z + dyz 
(doubly degenerate) and the band s+{dz2 ) for both spins. 
The band dz2 + (s) for the minority spin moves down with 
increasing U, and gradually becomes more d-like. At the 
same time, the dispersion and s character of the minor- 
ity spin s + (dz^) increases. In any case, two s channels 
are always open irrespective of U. The main difference 
between the band structure obtained with the LDA-I-SIC 
approach and the band structure for U—3 eV is that at 
the A edge, the minority spin band s + {dz^ ) almost do 
not move in case of LDA-I-SIC. For the other d-bands, the 
effect of SI is very similar to the results obtained with the 
LDA-I-C/ method using U of about 2-3 eV. 

Ni has one more electron than Co, and consequently, 
in the LDA-I-C/ approach, the dxy + dr^2_y2 (doubly de- 
generate) minority spin bands are now occupied and now 
shifted downward with increasing C/, in contrast to the 
situation in Co. As correlations increase within the d- 
shell, the majority spin channels with strictly d compo- 
nents are pushed downward. The partly unfilled minor- 
ity spin dxz + dyz (doubly degenerate) bands of Ni move 
slightly upward, while the corresponding majority spin 
bands move down. This is consistent with a small growth 
of the total magnetic moment Mt, (see Table lll|) . For 
the minority spin bands s + (dzs) and -|- (s) in Ni, one 
can see a similar change in dispersion to that already 
described for the Co wire. Also for this metal, two s 
channels are open at the Fermi level for all U studied. In 
the LDA-I-SIC method, the minority spin band + (s) 
moves a little bit downwards in comparison to the LSDA 
result, unlike the LDA-f C/ method. 

In Pd, the d^y + c?a;2_y2 bands move down with in- 
creased correlation within the d-shell, just as in Ni. Al- 
ready within LSDA, these orbitals are completely filled, 
making the magnetic behavior of this metal very differ- 
ent from the one of Co and Ni. A more relevant ef- 
fect of correlation in the d-shell is that the largely empty 
minority spin channel of the s -\- {d^2) band is pushed 
up, so that it becomes completely unoccupied already for 
C/=3 eV, whereas the corresponding majority spin chan- 
nel is pushed down. In the LDA+SIC band structure, 
the situation is similar. 

Furthermore, the dispersion of the minority spin chan- 
nel of the dz2 + (s) band increases so that this band be- 
comes increasingly s-like with larger U . For C/=3 eV or 
larger, only two bands cross the Fermi level, one minor- 
ity spin channel and one majority spin channel, both of 
s-|-d^2 character, with the s character dominating. Thus, 
also for this metal, the two s channels remain open in all 
cases. 

In contrast to LBA+U, the LDA-fSIC results yield 
a smaller spin polarization relative to LSDA within the 
d-shell, and similar magnetization within the s-shell as- 
sociated with the smaller total occupation of the s-shell. 

In Figure 01 we plot the s and the d components of 
the DOS for spin channels in the Pd wire. In both the 
LSDA and LDA+SIC calculations, the s-component for 
both spin-channels are present at the Fermi level. In the 
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LDA+U scheme with U =3 eV, the mmority spin at the 
Fermi level is dominated by the d-componcnt and the s- 
component almost vanish. For larger U , the d-component 
decreases and the s-component grows, as discussed ear- 
lier. 

Due to the aforementioned charge flow between the 
s- and d-shells, the magnetic moment in the s-shell de- 
creases because the s-occupation in the majority spin is 
constant and the s-occupation in the minority spin in- 
creases. While the d-occupation in the minority spin 
decreases and the magnetic moment within the d-shell 
increases. The above mechanism would keep the total 
magnetic moment constant, but the additional contribu- 
tion comes from the movement of d-bands from the Fermi 
level (opposite for the majority and minority spins). We 
stress that, the occupation within shells stems from the 
projected-DOS integrated over the whole energy range; 
therefore one should not conclude from Figure^ that a 
decrease of the projected-DOS at the Fermi level is corre- 
lated with a decrease of the occupation of the projected 
sheU. 

As for the effect of correlations for the stability, we 
calculated the equilibrium and breaking points within 
LDA+U with U=5 eV for all wires. The equilibrium 
distances (in A) for Co, Ni and Pd are 2.2, 2.11, 2.53 re- 
spectively, and the breaking points are 2.8, 2.63 and 2.96. 
The total-energy increase from the equilibrium point to 
the breaking point calculated within LDA/LDA-I-J7 with 
U=5 eV is (in eV) 0.751/0.666 for Co, 0.769/0.766 for 
Ni and 0.478/0.331 for Pd. The decrease of the energetic 
stability is mostly pronounced in Pd (about 30%) and 
the geometry change is the largest in the Co wire. This 
is somehow consistent with the band structure, because 
for Nickel only two bands disappear from the Fermi level 
when LDA-I-J7 with U=5 eV is applied while for Pd and 
Co wires three bands move out from the Fermi level. 

All above calculations were performed with the unit 
cells with one atom. We also calculated chains of seven 
atoms in the unit cell and with the U parameter in- 
creasing when going to the middle of a cell. This way, 
we wanted simulate somehow the "tips"; knowing that 
atoms at tip surface are less correlated than atoms in a 
chain and more correlated than atoms in bulk. 

In summary, for all three metals we find that both s 
channels (minority and majority spin) have finite frac- 
tional occupation for all values of U tested and also for 
the LDA-I-SIC calculations, which implies that at least 
the two conducting s channels are always present at the 
Fermi level. 



IV. DISCUSSION AND CONCLUSIONS 

In summary, we have found that inclusion of static 
repulsive electron correlations within LDA-I-J7, while al- 
tering noticeably the electronic structure of maximally 
stretched monatomic transition metal nanowires, does 
not affect radically the number and the nature of con- 



ducting channels. In particular, a decrease is found in 
the (poorly conducting) d channel number (Table II), 
whereas the calculated number of highly conducting s 
channels remains in all cases one per either spin. There- 
fore, one quantum Gq = 2e'^/h of s conductance will be 
approximately preserved after adding correlations at the 
static level of LDA+U. 

A certain reduction is predicted for the d channel con- 
ductance. At the same time, the nature of electronic 
states and their dispersion near Ep is also changed, 
which implies that tip-wire transmission should also be 
affected. The transmission changes are not addressed in 
this work, where tip-wire junctions are not explicitly in- 
cluded. However, because in LDA-I-J7 the filled states are 
lowered in energy and empty states are raised, the width 
of partly filled bands is somewhat larger, and that could 
in turn increase transmission. 

Nanowire magnetism (at zero temperature) is poorly 
affected by correlations in Co and Ni nanowires, whereas 
a clear increase is predicted in Pd. Antiferromagnetism 
(and the associated potentially insulating state) never 
prevails, at least for realistic U parameters and for inter- 
atomic distances below the estimated breaking point. 

It is problematic at this stage, where we have not in- 
cluded and treated at all the wire-tip junctions, to relate 
the present results to the experimentally reported ballis- 
tic conductance of Co, Ni and Pd nanocontacts. Never- 
theless we can qualitatively remark that the persistence 
of one s channel per spin leads to expect in all cases a 
conductance larger than one in units of Gq. While this is 
in agreement with break junction observations, where the 
ultimate conductance jump is around 1.3 in Co (Ref.lTol). 
1.4 in Ni (Ref.'T^, and 1.7 in Pd (Ref.HS), it does leave 
the fractional conductance steps G ~ 0.500 reported in 
STM experiments at room temperature completely unex- 
plained. More work, both experimental and theoretical 
will be needed to clarify this issue. 



Acknowledgments 

M.W. is grateful to Stefano de Gironcoli and Javier 
Junquera for discussion. We also wish to thank 
Ruben Weht for checking some of our results with 
the WIEN codaS, as well as D. Ugarte and C. Un- 
tiedt for discussions. Work in SISSA was sponsored by 
MIUR FIRB RBAU017S8R004, FIRB RBAU01LX5H, 
MIUR COFIN2003 and PRIN-COFIN2004, as weU as by 
INFM's "Iniziativa Trasversale Calcolo Parallelo". A.D. 
acknowledges financial support from VR (Swedish Sci- 
ence Foundation). 



APPENDIX A: SELF-CONSISTENT 
CALCULATION OF U 

We calculated U self-consistently for the wires. But 
also, we made a scan of the wire properties as a tunc- 
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tion of U, thus, using [/ as a free parameter. The main 
reason for this was that the calculated U comes out to 
be very large. When calculating U self-consistently, it is 
important to notice that the effect of U parameter on the 
band shift and on the other properties, in any implemen- 
tation of the LDA+[/ scheme, depends on the radius of 
the projector function, to, used to calculate the on-site 
occupation numbers, n^, (see Eq. ©). And U itself also 
depends on the projector radius, if it is calculated self- 
consistently. For smaller projector radius, calculated U 
is larger and the on-site occupation numbers are smaller. 
Therefore, the combined effect on the correction to the 
potential is more/less independent on the projectors. 

The self-consistent U cannot be directly compared with 
the Hubbard-?/, because the latter is defined directly by 
the spectroscopic properties, i.e. the electron affinities 
and the ionization potentials, as follows 

U = E{N + 1) + E{N -1) -2E{N), (Al) 

where E is the total energy of the system with A'' elec- 
trons in the considered shell. The above definition does 
not depend on the implementation. However, using the 
"exact" spectroscopic data for the U parameter could 
give a wrong picture, if such U was combined with the 
"not exact", i.e. projector dependent, occupation num- 
bers. Here we want to note that, for the spectroscopic 
definition of U we need to use atomic configuration which 
involves both s- and d-shelli^ Such calculations of U 
for atoms were performed for instance by Brandyopad- 
hyay and Sarma'^^ by means of the Hartree-Fock-Slater 
method. The phylosophy of standard LDA-|-?7 method 
applied in this work, however, is to correlate the elec- 
trons in the d-shell and not include the intershell interac- 
tions. Satisfying the above condition means that we can- 
not directly compare our U parameters with the atomic 
Hubbard-JJ, even if we abstract from the implementation 
and the combined effect of U and n„i- 

Focusing on the solid state methods, the Coulomb pa- 
rameter U can be determined from the linear response 
calculations with respect to a " suitable" perturbation in 
the occupation numbers (" suitable" means small but not 
too small, such that the method is sensitive to them) . 
The details of such an approach are given in Refs. jl7ll8l |. 
In this method, the U parameter is equal to 

U = ixo' - X-')h 

where index // denotes the diagonal element at a site /, 
which is the atom of our interest, and the noninteracting 
and the interacting response functions, i.e. xo and are 
defined as 

dnj Q dnj 

The quantity nj is the calculated local occupation num- 
ber of the considered polarization shell of the atom /. 
And a/ is the perturbation of the hamiltonian which en- 



TABLE III: Self-consistent U and the local occupations A^^j, 
as well as the Mulliken occupation numbers ]\jMuiuken 
the majority spin | and the minority spin J, within the d 
shell, calculated for the atoms in wires and for the separated 
atoms. 



parameters 


Ni 


Co 


Pd 


projector radius (a.u.) 


1.8 / 1.950 


2.073 


2.2 / 2.675 


Calculations for 


wires at the breaking point 


self-consistent U (eV) 


13.9 / 12.8 


12.7 


9.6 / 7.6 


U Nl (eV) 


56.0 / 53.5 


53.0 


34.8 / 31.8 


U A] (eV) 


39.1 / 37.5 


25.9 


30.8 / 28.5 


Nl 


4.03 / 4.18 


4.17 


3.62 / 4.18 


Nl 


2.81 / 2.93 


2.04 


3.21 / 3.75 


71 rt ,AIulliken 


4.91 / 4.91 


4.90 


4.89 / 4.89 


T\rl ,AIulliken 


3.53 / 3.54 


2.50 


4.39 / 4.43 


N'^ /A^'*^""'*'^" (%) 


82 / 85 


85 


74 / 85 


^j^i, Mulliken ^o^-j 


80 / 83 


82 


73 / 85 


Calculations for atoms (separation 


5 A in 


a chain) 


self-consistent U (eV) 


16.8 / 16.2 


15.6 


11.1 / 10.7 


U Nj (eV) 


70.4 / 70.3 


66.8 


40.4 / 44.5 


U Nl (eV) 


41.5 / 41.6 


29.3 


37.3 / 43.2 


Nl 


4.19 / 4.34 


4.28 


3.64 / 4.16 


Nl 


2.47 / 2.57 


1.88 


3.36 / 4.04 


T^r] ,A'Iulliken 
^^d 


5.00 / 5.00 


5.00 


4.98 / 4.95 


T\rl ,A'Iulliken 
^^d 


3.00 / 3.00 


2.19 


4.60 / 4.82 


^i/^]'*'"'"'"'" (%) 


84 / 87 


86 


73 / 84 


jyi ^j^l, Mulliken ^(y^-j 


82 / 86 


86 


73 / 84 



ters the total energy as follows: 

with the constrained occupation qj. The U parameters 
obtained this way are given in Table IHII for the Ni, Co 
and Pd wires and for the separated atoms. For Ni and 
Pd, we performed the calculations at two projectors radii. 

In Table Hill we collect occupations and calculated U 
parameters. As one can see, the value of obtained U 
parameter decreases with increasing radius of the pro- 
jector. The local occupation numbers within the d-shell, 
Nd, are smaller than the Mulliken ones, iV*^"'«'=™^ rpj^^ 
Mulliken occupation numbers were obtained by summing 
over contributions from all zeta functions in the basis set 
and they almost do not depend on the radius of the last 
(smallest) zeta- function (see Eq. (0J). While the local 
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occupations were summed over d-orbitals within the last 
zeta functions only, and these numbers depend strongly 
on the projector radius (see Eq. Q). The MuUiken oc- 
cupations sum to the total number of electrons in the 
unit cell, while the local occupations do not need to sum 
to total number of electrons within the c?-shell because 
they represent only the charge around nuclei and not in 
the whole unit cell. Our definitions given by Eq. ||2Jl are 
similar to the calculation of a charge closed within the 
atomic spheres in the LMTO approach (Ref. |lj|). The 
local occupations are smaller than the MuUiken ones, by 
typically 20% of Arj>^""ifc™_ 

In Table ITTll we give also the product UNd (for each 
spin) calculated at given projector radius. We assume 
that, the orbital occupations, n„j, are proportional to 
the averaged occupations, rid, and to the total occupa- 
tions, Nd, as well. Since the deviations of orbital occu- 
pations from the averaged ones, rim — nd, are typically 
one or two orders of magnitude smaller than the occupa- 
tion numbers rim and rid and Nd, the potential Vu (see 
Eq. ) is much smaller than the numbers UNd presented 
in Tab. Illll Nevertheless, the trend for the dependence 
of Vu and of UNd on the projector radius is the same, 
i.e. very weak. Thus, we may not bother too much what 
projector radius we chose, because the calculated prop- 
erties will not be much dependent on it, as long as it is 
in the reasonable range from 0.4 to 0.5 of the interatomic 
distance. 

The calculated values of the U parameters seem to be 
too large. But, we will show that, these numbers give 
an accurate estimation of the band shift calculated with 
the localized basis set code. We will use data for NiO (in 
the rhombohedral structure with the antiferromagnetic 
phase) where we know exactly the experimental band gap 
which is 4.3 eV fRef. ls^ . We calculated the fundamental 
gap in NiO with LSDA and LBA+U for U=2 eV and 
U^8 eV, and we obtained 0.22 eV, 0.69 eV, and 3.32 eV 
respectively, with the radius of the projector fixed at 2.5 
a.u. Fitting a parabola to above results, we see that, we 
should use the U parameter of about 9.5 eV to obtain the 
experimental band gap. Such value of the U parameter 
is already surprisingly high, not saying about values in 



Tab, mil if we compare them to other implementations of 
the LDA+C/ scheme. 

We give some reasons for using higher U than typi- 
cal values used in the plane wave approach. We rescale 
linearly the occupations for Ni obtained with the projec- 
tor radius 1.8 a.u. and 1.95 a.u. (see Tab. IIII|I to those 
(spin up and spin down) localized occupations which we 
would obtain using the radius of 2.5 a.u.; they would be 
4.85 and 3.31 for the majority spin and the minority spin 
respectively. Now, if we assume that the product UNd 
scales more/less linearly with the projector radius, we 
should take a value of about 48 for U TVj and of about 34 
for UN^ at the radius of 2.5 a.u. Dividing these numbers 
by the expected occupations at that projector radius, we 
will obtain the U parameters of about 9.9 eV and 10.2 eV 
for the majority spin and the minority spin respectively. 
Which gives in average about 10 eV. We know that the 
U parameter should be larger for the atom than for the 
bulk. Thus for the wire, a value of 10 eV is quite good in 
a comparison to its bulk value of 9.5 eV. This way, using 
the experimental data for a band gap in NiO and a fit 
for our numerical data obtained with LSDA and LDA-|-i7 
with U=2 and U =8 eV, and rescaling results for Ni wire 
to the expected result for a larger projector radius, we 
obtained two very close values of the U parameter for 
Ni. The result of above procedure is the proof that one 
could use the parameters presented in Table ITTll and that 
obtained this way band shifts in wires would be correct. 

On the other hand, the band shift is not the only prop- 
erty one would like to get from the LDA-I-J7 scheme. For 
instance, the correct orbital ordering in NiO could be ob- 
tained with a smaller value of U, about 6 eV. Therefore, 
for the sake of the balance between different properties 
which we want to describe, we take a bit smaller values 
of U than these which give the expected band shifts in 
bulk. 

At the end we conclude that, for the number of chan- 
nels crossing the Fermi level, it is not very important how 
large U we take, because for U about 3'^4 eV the bands 
which are supposed to move already have done it. 
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